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ABSTRACT 

We present the results of the decade-long M31 observation from the Wendel- 
stein Calar Alto Pixellensing Project (WeCAPP). WeCAPP has monitored M31 
from 1997 till 2008 in both R- and J-hlters, thus provides the longest baseline of 
all M31 microlensing surveys. The data are analyzed with the difference imaging 
analysis, which is most suitable to study variability in crowded stellar helds. We 
extracted light curves based on each pixel, and devised selection criteria that are 
optimized to identify microlensing events. This leads to 10 new events, and sums 
up to a total of 12 microlensing events from WeCAPP, for which we derive their 
timescales, flux excesses, and colors from their light curves. The color of the 
lensed stars fall between {R — I) = 0.56 to 1.36, with a median of 1.0 mag, in 
agreement with our expectation that the sources are most likely bright, red stars 
at post main-sequence stage. The event FWHM timescales range from 0.5 to 14 
days, with a median of 3 days, in good agreement with predictions based on the 
model of Riffeser et al. (2006). 

Subject headings: dark matter — gravitational lensing — galaxies: halos — galaxies: 
individual (M31, NGC 224) — Galaxy: halo — galaxies: luminosity function, mass 
function 
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Introduction 


Almost four decades after revealing evidences for dark matter in galaxies flRubin et al 


1980), its nature is still unknown. Dark matter can be smoothly distributed, e.g. the 


weakly interacting massive particles (WIMPs), or in compact form. Since dar k matter do es 


not em it light, the best way to study it is through gravitational interaction. 


Paczvnski 


(119861) was the hrst to conceive the idea of gravitational microlensing as a method to 


detect massive compact halo dark matter (MACHO). Based on his calculation, the optical 
depthlj, i.e. at any time the probability of a source to be closer along the line of sight to 
a foreground lens than the lens’ Einstein angle therefore to be magnihed by more than 
1.34, is of order 10“® towards Magellanic Clouds. This has motivated seve ral campaigns to 


search for microlensing towards dense stellar helds, see e.g. t he review by 


lie first micro 


(lAubourg et al. 


lJ 


ensing events were reported by the MA CHO (lAlcock et al. 


1993 ). and OGLE f Udalski et al. 


Vloniez 


1993) teams, with MO A ( Muraki et al. 


(120 inh . 


99311. EROS 


19991 ) joining later on. In their 5.7-year survey data, the MA CHO team has an nounced 


microlensi ng events towa rds the Large Magellanic Cloud flAlcock et al. 


later analysis (jBennett 


20001) . while a 


20051 ) showed that 1 of them is a true variable star and 2 of them 


are likely to be variables from a simple likelihood analysis. This leaves 10 of them to 


be plausible microlensing ev ents, anc 


masses between 0.1-1 Mq ([Bennett 


gives a MACHO halo fraction of 16% for MACHO 


2 OO 5 I) . The EROS-1 and EROS-2 surveys resulted 


in a MACHO halo fraction of less than 8% for MACHO mass ~ 0.4 Mq, and ruled out 
MACHOs with masses between 0.6x10“^ and 15 Mq as a major component of the Milky 


Way halo fjTisserand et al.l 


20071) . Analyzing the OGLE phase H and HI data, Wyrzykowski 


^We note that this valu e depends primarily on the lens population characteristics. 


In flAlcock et al. 


20001) they found 13 microlensing events with tight criteria; they also 


presented 17 events with loose criteria, which contain more low S/N events. 
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et al. (2010,2011a,2011b) concluded that microlensing events towards the Large and Small 
Magellanic Cloud can be reconciled with self-lensing by stars alone , i.e. w ithout requiring 
compact halo objects in the Milky Way halo. Later on. 


Besla et al 


fl2013l) studied the tidal 


streams between the Magellanic Clouds and used theoretical modeling to show that the 
microlensing signal can be r eproduced by the stars in the stream, though their modeling 


requires further verihcation. 


Besla et ah 


(120131) outlined several observational tests to verify 


their theoretical modeling, e.g. the sources of the microlensing events are low-metallicity 
SMC stars, the sources have high velocities relative to the LMC disk stars, and the presence 
of a very faint stellar counterpart to the Magellanic Stream and Bridge (surface brightness 
> 34 mag/arcsec^ in V-band). 


In contrast to this. 


Calchi Novati et ah 


(120131) recently re-analyzed the OGLE and 


EROS data and found that some of the OGLE events can be attributed to halo lensing. 
The inconclusive results can be partially attributed to the fact that, by monitoring the 
Magellanic Clouds we can only sample a small fraction of the Milky Way halo. In order to 
have a census of the Milky Way halo, we need targets that are distributed along different 
lines-of-sights through the Milky Way halo. Other than the Milky Way, we could monitor a 
nearby spiral galaxy, for which we can have a complete view of its dark halo. 


An alternative dense ste 
proposed by several authors i 


lar target for a microlensing search for MAC H Os is M31, as 


Crotts 

1992; 

Baillon et al. 

1993; 

Jetzer 

1994) 


19941) . The advantage 


of M31 is twofold. First, we can have multiple lines-of-sight through a dark halo towards 
M31, in contrast to the single line-of-sight towards the Magellanic Clouds due to our location 
in the Milky Way. The second advantage is that one does not only probe the Milky Way 
halo, b ut also that of M31. The stru cture of M31 is well known: with an inclination angle 


of 77° (jWalterbos fc Kennicntt 


19871) . we expect to see an asymmetry of the microlensing 


event rate between the near- and far-side of M31 disk. Such an asymmetry can however 
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also originate from extinction along different sight-lines caused by dust. This has indeed 
been observed i n the density distribution of variables from the POINT-AGAPE survey 


flAn et ah 


2004 1 and from the We CAP P f Fliri et al 


dust effect, we have studied the dust properties of M31 flMontalto et ah 


2006). To be able to account for this 


20091) and derived 


an extinction map across the disk of M31, which can be used for quantifying detection 
efficiencies. Micr olensing in M31 differs from the Magellanic Clouds, because at the distance 


of M31 (770 kpc. 


Freedman fc Madore 


19901) . most of the sources for possible microlensing 


events are not resolved, and each pixel of a CCD image contains up to hundreds of stars. 
Instead of monitoring individual resolved stars, one has to monitor pixel light curves and 
their variations towards dense stellar helds. The sources of microlensing events are usually 
not resolved before and after the high magnihcation phases and their “baseline fluxes” are 
thus unknown. In most cases the magnihed source (at least at maximum magnihcation) 
appears as a resolved object. However, in order to achieve exquisite photometry in such 
crowded helds, we thus perform diherence image analysis and search for microlensing events 
at the position of each individual pixel. In order to extract light curves at the position 
of each individual pixel, we use a PSF constructed from resolved sources in the image to 
perform PSF photometry and obtain pixel light curves. 


The theoretical aspect of p ixel lensing ha s been laid down by 


Riheser et al. 


Gouldl (119961) under small 


(120061) reformulated this theory in a more 


impact parameter assumption, 
general context. 

The hrst microlensing events towards M31 were re ported by the VATT/Columbia 


micro lensing survey. They put the idea of 


CrottsI (119921) into practice (ITomanev fc Crotts 


19961) and presented 6 microlensing events discovered by the joint observations of the 


Vatican Advanc ed Technology Telescope (VATT) and the KPNO 4m telescope taken during 


1994 and 1995 ( Crotts fc Tomanev 


19961) . Their observations continued from 1997 to 
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1999, with the VATT and the 1.3m telescope at the MDM observatory. With additional 


data from the Isaac New ton Telescope, 


out of their 3 year data 


Uglesich et ah 


hey 


presented 4 probable 


2004 ). At the same time. 


microlensing events 


Ansari et ah 


1119971) 


also launched the Andromeda Gravitational Amplihcation Pixel Experiment (AGAPE); 
they observed M31 with the 2m telescope Bernard Lyot (TBL) in the French Pyrenees 
i n 1994 and 1995, which led to the discovery of one bright, short microlensing event 


flAnsari et al, 


19991) . Following AGAPE, the Pixel-lensing Observations with the Isaac 


Newton Telescope-Andromeda Galaxy Amplihed Pixels Experiment (POINT-AGAPE) 


Auriere et al. 

(2001 

), accompanied bv three more in 

Paulin-Henriksson et al. 

(2002, 

2003) 


another 3 were reported by 


Galchi Novati et al. 


(120031) with additional data from the 1.3m 


telescope at the MDM observatory in 1998-1999. The full POINT-AGAPE data were 


Belokurov et al. 

2005 

), 6( 

Galchi Novati et al. 

2005 

), and 10 events ( 

Tsaoras et al. 

2010) 


were reported by the different nodes, respectively. 


Using the same INT data, the Microle nsing Explora 


(MEGA) survey has presented 14 events fide Jong et al. 


ion o f the Galaxy and Andromeda 


20061) . At the same time, the 


Nainital Microlensing Survey employed the 1.04m Sampurnanand Telescope in India to 


observe from September 19 98 until Fe 


event from the 4-year data fjJoshi et al 


jruary 2002. They have extracted 1 microlensing 


20051) . More recently, the Pixel Lensing Andromeda 


collabor ation (PLAN) carrier 


in Italy fiGalchi Novati et al 


(iGalchi Novati et al. 


out observations using the 1.5m Loiano telescope located 


20071) and reported 2 events in their data collected in 2007 


20091) . They further incorporated observations fro m the 2m Himalayan 


Chandra Telescope (HCT) taken in 2010 and reported another event (jCalchi Novati et al. 


20141 ). In the mean time, the Pan-STARRS 1 collaboration conducted a high-cadence. 


long-term Andromeda monitoring campaign (PAndromeda) utilizing its wide-held (~ 7 
































































deg^) camera. Based on its first year data, iLee et al.l fj2012l) have reported 6 events in the 
central 40”x40” area, with the promises to detect more events taking advantage of the 
larger survey area and higher cadence from Pan-STARRS 1. 

It is worth to note that previous M31 microlensing event identih cations may suffer from 


contaminations by variables. For example. 


Crotts fc TomanevI (119961) suspected that part of 


their events are contaminated by long-period red supergiant variables. Despite the efforts of 
vari ous campaigns, the M ACHO fraction at the mass range of O.I-IMq is still under debate 


see 


Calchi Novati 


2012 


, for a detailed discussion). For 


example, the POIF 


collaboration has reporte d evidence for 


the MEGA collaboration fide Jong et ah 


MACH O signal f Calchi Novati et al 


T-AGAPE 


20051) . while 


20061) on the contrary concluded that their events 


can be fully explained by self-lensing. 


Due to the small number of reported events, the origin of M31 microlensing remains 
an open issue. In this study, we aim to increase the number of microlensing detections and 
suppress contaminations by variables with long-term observations of the M31 bulge. This 
paper is structured as follows. In section H] we present the observations of our long-term 
survey. Our data reduction is outlined in section [3l followed by the event detection in 
section 01 The analysis of these events are shown in section |5l A discussion of our events, 
as well as results from previous M31 microlensing surveys are presented in section [61 with a 
summary and prospects in section [71 


2. Observations 


WeCAPP continuously monitored M31 from August 1997 until March 2008 using the 


Wendelstein 0.8m telescope flRiffeser et al 


200ll) . The data were initially taken with a TEK 


CCD with 1,024x1,024 pixels with a £eld-of-view of 8.3x8.3 arcmin^ pointing at the bulge 





























of M31, optimally o n a d aily basis in both R- and /-filters. Following the suggestions of 


Tomanev fc CrottsI f 19961) and 


Han fc Gouldl f 19961) . we pointed to the far side of the M31 


disk (FI in Fig. [T]), where the halo lensing probability is maximized. From June 1999 to 
December 2002 we collected additional data using the 1.23 m (17’.2x17’.2 FOV) telescope 
at Calar Alto Observatory in Spain to increase the time sampling. This provided a FOV 
which is four times the Wendelstein FOV, and enabled us to survey the major part of the 
M31 bulge. After 2002, we used the Wendelstein telescope solely to mosaic the full Calar 
Alto £eld-of-view with four pointings, as indicated in Fig. [1] 


The amount of observations taken in the four pointings differs signihcantly during the 
11 seasons. Fig. [2] shows a histogram of the number of observed nights by WeCAPP. The 
most complete seasons are 2000/2001 and 2001/2002 with joint observations from both 
Wendelstein and Calar Alto (see also Tabled]). 


To have an overview of the observing cadence, we show the daily sampling in Fig. [S] 
Thanks to the joint observations at Wendelstein and Calar Alto, we achieved an average 
time coverage for FI in R of 42% during the 2000/2001 season (peaking in August 2000 
with 90% on JD ~ 2451770) and an average time coverage of 55% during the 2001/2002 
season, reaching more than 93% in 3 months (July and October 2001, and January 2002, 
around JD ~ 2452110, 2452200 and 2452290, respectively). During the 11 seasons, we have 
obtained a total sampling efficiencjo of 14.9% in R and 11.5% in /, with 11.3% for R and I 
combined for all helds (F1-F4), which means that in 11.3% of the nights we have both R 
and I observations. 


However, not all images have the same quality. Rather than quantifying the fraction 


^The sampling efficiency is a fraction with respect to the overall baseline length, i.e. 
including periods where no observations were scheduled. 
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Fig. 1.— M31 composite image {V-, R-, and /-band) of the observed fields FI to F4, taken 
at Calar Alto Observatory dnring the 2000/2001 season. The black lines mark the positions 
of helds FI to F4. The Calar Alto camera covers all fonr helds simnltaneously. Field FI was 
observed during all 11 campaigns from September 1997 until March 2008. 

of nights we have observed through the 11 seasons we would like to have (as a function of 
location in M31) the fraction of nights where the noise is below a certain threshold. For 
this reason, we empirically chose a noise limit of the minimum S/N of 8.9 for our faintest 
event, i.e. 0.73 x 10“^ Jy/8.9, which corresponds to the 8.9(T detection criterion we present 
in section 4. For pixels with noise levels above this value, the lensing signal would mix 
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R-band I-band 
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JD - 2450000 JD - 2450000 

Fig. 2.— Histogram of the number of analyzed nights. The 4 different helds FI, F2, F3, F4 
are colored blue red, yellow and green, respectively. Periods marked in gray show the 61 days 
(from 1®* April to 3H* May) during which M31 can hardly be observed (see also Table. [I]). 
Our most complete seasons (2000/2001 and 2001/2002) are marked with the texts “00-01” 
and “01-02”. 

with the high noises and could not be detected. Hence these pixels cannot be used for the 
detection. The sampling thus depends on the x and y position of the pixel, as well as the 
observation time t, which we denote as < S{x,y,t)>. In Fig. IH we show the area having 
a noise smaller than our noise limit for every observed night. By averaging over time t, we 
will get the positional dependence as shown in Fig. [5l 


Data reduction 


We process the data using our customized pipeline MUPIPE flGossl fc Riffeserl 120021 ). 


where standard reduction processes - such as bias subtraction, treatment of bad pixels, 
flat-helding, cosmic ray removing - are performed with per pixel error propagation. To 
identify unresolv e d var iables, we employ the difference imaging analysis (DIA) proposed by 
Alard fc Luptonl (119981) . which enables us to detect variables with amplitudes at the photon 
noise level and to measure their flux excesses relative to high signal-to-noise reference 
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R-band 



JD - 2450000 


Fig. 3.— Daily sampling of the 4 different fields FI, F2, F3, F4 are colored in blue, red, 
yellow and green, respectively. Periods marked in gray show the 61 days (from 1®* April to 
31^* May) during which M31 can hardly be observed. See also Table. [T] The two seasons 
with the highest sampling are 2000/2001 and 2001/2002, i.e. those seasons where we could 
combine Wendelstein with Calar Alto data. 


images. 


After the difference imaging, we perform PSF photometry on each pixel in the following 
manner. We hrst extract the PSF prohle from several isolated, bright and unsaturated 
reference stars. Then we £t this PSF to all pixels to generate light curves of varying sources 
identified in the difference images. The flux of the source is estimated by integrating the 
count rates over the area of the PSF. 


(12003 


le res ults from a subset of data of this project have been pres ented in iRiffeser et al 


20081) and partially contributed to 


Calchi Novati et al. 


(2010). In addition to the 


original microlensing targets, the high cadence observations also yielded a sample of more 
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JD - 2450000 


Fig. 4.— Sampling in t averaged over x and y < S{x,y,t) > 


xy 


than 20,000 variabl 

es in the 

mlge 

IPietsch et ah 

2007; 

Lee et ah 

2012) 


of M31 f Fliri et al.l 12006 1 and 91 candidate novae 


Event detection 


Paczynskil fjl9861) provided an analytic formula to describe the amplification of a 
microlensing event: 


A{t) = 


u^ + 2 


u = 






( 1 ) 


u\/u^ + 4’ 

where to is the time of maximum amplihcation, tg is the Einstein ring crossing time, Uq is 
the impact parameter in units of the Einstein ring radius. The light curve (or the measured 
flux as a function of time) of the microlensing event can thus be expressed as: 


F{t) = F^[A{t)-l]+B, 


(2) 


where Fq is the un-lensed flux and B is the blending within the PSF. 
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Fig. 5.— Sampling in {x,y) averaged over time < S{x,y,t) >t- Color levels range from 0% 
to 10%; the two contour levels show 1% and 5%, respectively. 58 arcmin^ are well sampled 
by more than 5% of the 11 year survey. 220 arcmin^ are well sampled by more than 1%. 

This conventional microlensing light curve formula is highly degenerate in and uq 
for microlensing events towards M31, because one can not resolve the source anymore. For 
M31 microlensing, the only observables are the flux excess and the event time-scale 
^FWHM- Gould (1996) has derived the pixel-lensing light curve formula in the context of high 
magnihcation. Riffeser et ah (2006) further revised the microlensing light curve formula 
by Gould (1996), such that microlensing light curves with moderate magnihcations can be 
well-described as well. The formula of pixel-lensing light curve from Riffeser et ah (2006) is 
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expressed as: 


Apit) 


eff 


m-tof 

''FWHM 


1 - 1/2 


+ 1 


B 


(3) 


where Fgg is the effective ffux, which for high magnifications is approximated by the flux 
excess F^s = Fq/uq ~ Ap. Since Paczynski (1986) was the first to present the simple 
analytical form of microlensing light curve as equation (|2]), we refer to equation ([2]) as 
Paczynski-fit throughout this paper. Since Gould (1996) was the first to introduce the 
approximated form of microlensing light curve as equation ([3]), we refer to equation ([3]) as 
Gould-fit throughout this paper. 


We use equation (E]) to identify microlensing events. The microlensing event detection 
is performed on the light curve of each pixel (based on the aforementioned PSF photometry) 
with several successive criteria. Gompared to Riffeser et ah (2003), we introduce additional 
criteria to avoid human interaction during the selection process. We now describe our 
criteria; an overview of the number of pixel light curves passing each criterion is listed in 
Table |2l 


• Griterion I is applied to exclude pixel light curves which have too few data points to 
make them worth analyzing. Hence we exclude light curves with less than 50 data 
points in either R- or I-filters. This leaves us with 3,872,240 light curves. 


• Griterion II identifies varying sources warrant for further analysis. We preselect 
variable pixel light curves with at least three consecutive 3a outlier^ in the R-filter. 


'^The difference between the individual flux measurement and the constant part/offset of 
the light curve is three times larger than the error of this individual flux measurement. The 
constant part/offset of the light curve is derived in two steps: after removing all data points 
with very high errors (larger than 5 times the median of all errors) we determine the constant 
part/offset of the light curve by iteratively fitting a line and clipping away all points which 
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• Criterion III is designed to find microlensing light curves with good We select 
pixel light curves that are well described by microlensing light curves. We use the 
microlensing light curve in equation ([3]) with some modifications: 


A pit) ^ Feflf 


(4 + 0 . 5 )^ 


- 1/2 

+ C 


(4) 


where we set tgt = i/tFWHM — 0.5, hence the full width half-maximum event timescale 
is always greater than or equal to 0.5 day. We use tat instead of fpwHM to avoid 
un-physically small fpwHM due to limited time resolution (up to half day) of our 
observation. The term c takes into account the shift of the baseline flux by a constant 
in cases for which there is a variable source spatially close to the microlensing event 
and data of these variable phases enter the reference image. 


We then filter out light curves with a Xk = \/xf (where Xk total divided by 

the degrees of freedom k) larger than 1.5 in the R band and larger than 2.1 in the I 
band. This is less strict than the value we used in (Riffeser et al. 2003). The allowed 
Xk in I is slightly higher than R, because the noise level is increased by unaccounted 
systematics such as the detector fringing and nearby variable stars in the I band. 


• Criterion IV evaluates whether a high flux excess relative to the baseline is more 
likely caused by random noise or due to a true microlensing signal. We consider 
the S/N of such a high flux excess measurement, SN* = , where the flux offset 

is Acon,i = Vi — c, Hi is the ith flux measurement at time ti and cTj is the error of 
the flux measurement. We take into account the probability of such a high flux 
excess measurement being close to the maximum of a microlensing light curve fit 
oc exp where the flux offset is = i/i — ApiU), and ApiU) is the 


are more than 2 sigma/errors above. 
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model flux at time t* according to Equation (|1]). We also take into account the 
probability of such a high flux excess measurement being at the constant baseline 



We then combine these three factors and define the S/N 


Pcon,i oc exp 
probability (SNP): 


SNPi = SNj X X (1 - Pcon,j) 


(5) 


This ensures that light curves that have high S/N outliers, that are outside the time 
interval of the microlensing event, get a lower weight. 

To avoid multiple detections of the same microlensing signal in neighboring pixels, we 
only use the pixel with good PSF detectionj^ to evaluate SNP,, i.e. where the PSF fit 
has a minimum in with respect to neighboring pixels. 

We empirically require one data point in the light curve to have a SNPj larger than 
8.9 in the R band to efficiently reject faint variable sources. 

• Criterion V quantifies the temporal correlations of the model-data mismatch (best-fit 
microlensing light curve vs. measured light curve) and enables us to reject intrinsic 
variable sources more efficiently. 

We combine the probabilities for positive p// or negative flux offsets pj" from the 
best-fit microlensing light curve, and assign positive values for consecutive data points 

®To find the good PSF detections we have carried out the following steps: 1. We subsam¬ 
ple each pixel by a factor of 5; 2. We determine the flux at the position of each sub-pixel by 
fitting a PSF; 3. We compare the of fh® PSF with the adjacent 8 pixels. If the pixel has 
the lowest among the adjacent 8 pixels, we consider it as a good PSF detection. After 
the detection we refit the position of the microlensing events and determine their positions 
at sub-pixel level. 




that are always on the same side (either above or below) of the best-£t model light 
curve. For consecutive data points alternating along the best-£t model light curve, 
we assign negative values. We then dehne the energy of a potential microlensing light 


curve as 


n—1 



( 6 ) 


where n is the number of data points and the probabilities are dehned as 



= 0.5 1 + erf 


and Pi =l-p+. 


For a random process the distribution of F' is a Gaussian with an expectation value 
of zero and a standard deviation of one. 


We derive the energy using the n = 20 closest data points to to in the R-band. We 


hnd empirically that a value between -1.5 and 1.5 efficiently rejects periodic variables 
and allows to skip the previously used by-eye detection in Riffeser et al. (2003). 

Criterion VI biters out long-periodic variables. Since our light curves span a baseline 
of 11 years, contaminations from the long-periodic variables are less severe than any 
of the previous campaigns. We inspected the candidate light curves and found that 
false detection from systematics and moving objects are having timescales longer 
than 1000 days. Hence we are able to empirically increase the upper ipwEM-hmit for 
micrlensing searches to 1000 days (compared to, e.g. trwHM < 20 days in Riffeser 
et la. 2003). With this criterion we also filter out objects with proper motion. A 
moving object that passes through the field with constant angular velocity will cause 
variabilities at some (stationary) pixels that can mimic microlensing signals. These 
proper motion objects can of course be excluded by inspecting postage stamps, but 
we would like to have a completely automatic selection here. 


Criterion VII rejects light curves which look like microlensing events from their overall 
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light curves but which are not well sampled close to the light curve maximum, hence 
could be variable sources. The contaminations are mostly from novae, e.g. not well 
sampled in their rising parts. 

We dehne the sampling quality for the falling and rising parts of each light curve 
within (to —15 days, to) and {to, to +15 days). The contribution of a single data point is 
then calculated by integrating the model light curve within (fj — 0.5 days, ti + 0.5 days). 
As sampling criteria we require a total sampling of the area under the light curve of 
at least 18% on the rising part of the light curve and of at least 8% on the falling side 
in the R and I band. 

In the 4*^ column of Table |2] we show how applying each criterion (II, III, IV, V, VI, 
VII) to the 3,872,240 light curves (criterion I) reduces the hltered light curves. For example, 
152,753/3,872,240~3.9% light curves that pass criterion III. Less efficient criteria show a 
high percentage of remaining light curves while efficient criteria hlter out more light curves. 
Therefore the criterion IV is our most efficient criterion. In the 3rd column we show how 
applying each criterion (II, III, IV, V, VI, VII) to the 3,872,240 light curves (criterion I) 
reduces the hltered light curves. 

Fig. [6] shows the to distribution of all pixel light curves under analysis. The colored 
histograms show how the number of light curves is reduced after each criterion. The total 
numbers correspond to those given in Table [21 

Because of large dome seeing and an inappropriate auto-guiding system, photometric 
errors are largest during the hrst season (1997/98). During the second season (1998/99) 
we were able to decrease the FWHM of the PSF by a factor of two, and therefore the 
photometric scatter also became smaller. 

The vast majority of our microlensing events are detected between 2000/2001 and 
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Julian date - 2450000 


Fig. 6.— Histogram of the number of light curves passing different detection criteria. Our 
most complete seasons (2000/2001 and 2001/2002) are marked with “00-01” and “01-02”. 
Yellow: x-limit (HI). Black: {S/N) constraint and good PSF (IV) at light curve maximum. 
Green: Energy criterion for the 20 closest data points to the maximum (V). Blue: timescale 
constraint (VI). Red: sampling criterion for the closest data points to the maximum (VII). 
In the end we are left with 12 microlensing events detected (see Table [2] for details). Eight 
of them took place during the two best observing seasons 2000/2001 and 2001/2002. 


2001/2002 seasons, where we have employed the Calar Alto telescope and the cadence is 
high. This implies that deeper images with larger telescope along with densely sampled 
observations are crucial to detect microlensing events. Our observations in other seasons 
are pivotal as well, as they serve the purpose to rule out contamination from variables. 


Fig. [7] shows the properties of the 719,628 light curves passing criterion II. Since the 
major contaminations to microlensi ng detections are variables and novae, we over-plot the 


23,001 variable sources published in 


Fliri et al. 


fj2006f) in magenta and the 91 novae in Lee 


et al. 2012 in blue. In Fig. [7] it is clear that criterion III and IV are efficient to hlter out 
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these two major contamination sources. 



inax(XR,%T-0.6) 


Fig. 7.— Criterion III and IV reduce the light curves from 719,628 (green) to 2247, as 
shown in the box on the upper left side. The light curves which do no pass criterion VI 
are marked in yello w. We show the 23,001 variables published in the WeCAPP variable 


star catalogue from 


Fliri et ah 


(120061) in magenta, as well as the 91 novae from Lee et ah 


(2012) in blue. The hnal 12 microlensing events are shown in red. We note that the decline 
parts of nova light curves resemble microlensing events and often make short timescale novae 
the major contaminations to microlensing detections. In this regard, we show that criterion 
III and IV can efficiently discriminate microlensing events from novae, hence remove their 
cont aminat ions. 


We note that not all of the criteria are independent. To disentangle the issue of criteria 
overlapping with each other, we perform a test with a subset of the criteria. The results are 
shown in Table [31 

Grit II and IV overlap already by dehnition. All detections passing Grit IV are a subset 
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of the detections passing II, because we fit only light curves which have 3 times 3-sigma 
outliers (Grit II). As can be seen in row (3) only 2.6% of the Grit II light curves pass the 
Grit IV, therefore Grit IV is much stricter than Grit II. A comaprision between row (7) and 
(10) shows that Grit IV is able to reduce the preselected light curves from all other criteria 
by a factor 3400. 

The comparision between Grit III and V shows that Grit V is weaker than III, because 
it reduces the light curves only by 56% compared to 21% from Grit III. Two-thirds of the 
light curves passing Grit II overlap with the light curves passing Grit V (14.8% of 21.2%), 
this shows a quite strong overlap. Also row (8) shows that Grit V has only a small impact 
to the detection as it reduces the number of detection only from 15 to 12. 

In short, we are aware of the overlapping among criteria. However, in the cases of Grit 
II and IV, we need Grit II prior to Grit IV to preselect the light curves of interest, so that 
we can analyze the full data-set in a reasonable amount of time. Nevertheless overlapping 
criteria should not be an issue if the efficiency study is running through the same detection 
criteria. 

Note that in Griterion VI we deliberately set the upper limit of fpwHM to be 1000 days. 
The fact that our hnal set comprises only events with tpwHM smaller than 15 days, can be 
hint that no rather long microlensing events in M31 exists. If we would have set the limit 
to 20 days, we could not address this speculation. 

This tpwHM limit is empirically picked by looking into the fpwHM distributions of the 
2247 detections fulhlling passing criterion IV, as shown in Fig. [HI 

We note that the long timescale hts may arise from slowly moving objects. To test 
for this we inspect the positions of all detections passing criterion IV with tpwHM > 1000 
days. We hnd that they are homogeneously distributed (they appear in some spikes of 
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Fig. 8.— tpwHM distrubtions of the 2247 detections passing criterion IV. The blue vertical 
line indicates the 1000 days upper limit in our selection criterion. 

bright stars) except a small region close to a bright high proper motion object. This object 
influences the difference imaging kernel and therefore mimicking a proper motion of all 
other sources. 

To demonstrate that 1000-d is a reasonable limit for fpwHM, we plot the fpwHM and 
distributions of the events in Fig. [9l The yellow points are all detections passing criterion 
IV and with fpwHM > lOOOd. The green points are all detections passing criterion IV and 
with fpwHM < lOOOd. The magenta points are variables from Fliri et ah (2006). The 
blue points are novae from Lee et ah (2012). The red points with black dot are the 12 
microlensing events presented in this paper. 
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X 



0.9 1.0 


1.1 1.2 1.3 1.4 1.5 

max(XR,Xi-0-6) 


1.6 


Fig. 9.— The events passing criterion IV, as well as our selection criteria of tpwHM and x^- 
The yellow points are all detections passing criterion IV and with tpwHM > lOOOd. The green 
points are all detections passing criterion IV and with tpwHM < lOOOd. The magenta points 
are variables from Fliri et ah (2006). The blue points are novae from Lee et ah (2012). The 
red points with black dot are the 12 microlensing events presented in this paper. 

5. WeCAPP M31 lensing events 

In this section we present the 12 microlensing events found in the WeCAPP data. 
Table 0] gives an overview on the properties of these 12 events. In this table we focus on 
two observables fpwHM and Ap. We fit the light curves with equation [2l 

As shown in Table IH the Paczynski-fit provides additional information on the Einstein 
timescale tp and the impact parameter uq. 


We note that the parameters tpwHM and Ap are degenerate, which often leads to 
overestimated errors in these two parameters for very short time-scale events (tpwHM below 
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2 day). Irrespective of this degeneracy, we can still determine these two parameters with 
reasonable nncertainty. To demonstrate the degeneracy, and onr conhdence in determining 
these two parameters, we show the 1, 2, and 3-a contonr plot of fpwHM vs. Ap in the 
appendix. 


The light cnrves of the 12 WeCAPP events are shown in Figs [TUlfT^ To illustrate the 
microlensing nature we also present postage stamps from the difference images close to the 
light curve maximum. These images help to rule out artifacts overseen by the pipeline, like 
hot pixels or cosmic rays as origin for the events. Indeed, as the inspection of the postage 
stamps show, none of our events is due to such an overseen artifact. In some cases (W2, 
W5, W7) the postage star nps also show positi ons of nearby variables. Because of improved 


photometric methods (see 


Riffeser et al 


differ from the published ones in 


2006h. the li ght curves for W1 and W2 slightly 


Riffeser et al. 


(120031 ). but agree within the error bars. 


6. Discussion 

In this section we summarize M31 microlensing events from previous surveys we are 
aware of and we compare our 12 events with previous studies. 

Table |5] lists the project name, number of detected events, their label and color as used 
in Figures ITdl and ITdl and the references for the events. In the 2"'^ column we report the 
number of events that were solely detected by the corresponding survey. In the 3’’'^ column 
we list the number of events that were detected by more than one groups. For example, 
Paulin-Henriksson et al. (2003) published 4 events from the POINT-AGAPE project, where 
2 of them were only detected by POINT-AGAPE, and 2 of them were also detected by 
other projects. Altogether 56 different M31 microlensing events have been published to 
date (including those presented in this work). Their fpwHM, Ap distributions are shown 
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in Fig. [131 In Fig- HHwe show the positions of microlensing events in a 50x70 arcmin^ held 
from the 2"''^ Palomar Sky Survey. 


The hrst fact to note in the middle panel of Fig. [13] is that WeCAPP (red), 
POINT-AGAPE (blue), and PLAN (yellow) surveys detected very short timescale events 
and hardly any event s with ^fwhm > 15 days. In contrast, the hrst events ever announced 


by VATT/Columbia f Crotts fc Tomanev 


1996 


, marked in violet color in Fig. [T3ll have very 


l ong t imescales (5 of 6 events have tpwHM > 65 days). In the same paper 


Crotts fc Tomanev 


(119961 ) suspected that part of their events may be caused by long period variables. The 


MEGA (green) and Belokurov L2 (cyan) events show an almost hat timescale distribution. 
Five of the MEGA events are (in de-projection) located in the disk, almost at the same 
distance with respect to the M31 center. This makes them more suspicious, since they 
could be intrinsic varying disk stars that exist in a particular evolutionary time scale, that 
is overrepresented at this particular radius within the disk. 


Predicting event rates and their characteristics from the theoretical studies, Riheser 
et ah (2008) from their M31 microlensing model show the expected event rate of long 
timescale events (tpwHM > 30 days) is an order of magnitude smaller than the event rate 
of short timescale events (tpwHM < 10 days), see e.g. their Fig. 12. The fact that the vast 
majority of the events presented in this paper are short timescale events is therefore in good 
agreement with expectations. In contrast, VATT/Columbin and MEGA surveys show an 
almost flat timescale distribution. Unless surveys analyzed by MEGA and Belokurov et ah 
(2005) have detection efficiencies that strongly suppress events with timescales < 10 days, 
they should have detected more short timescale microlensing. An alternative and more 
plausible explanation is that a large fraction of these long events are indeed long period 
variables, but were not ruled out from the microlensing candidate light curves due to the 
relative short time-span of the survey. 
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The VATT/Columbia team uses their observations in the 1994-1995 season to rule out 


obj ec ts of masses in the 


range o 


M31 f^Crotts fc Tomanev 


1996|). 


0.003-0.08Mf7| as the primary constituents of the mass of 


Uglesich et af 


(120041) hnd 4 probable microlensing events 


with data collected between 1997 and 1999. They conclude that 29^^°% of the halo masses 
are composed by MACHOs assuming a nearly singular isothermal sphere model. They also 
provide a poorly constrained lensing component mass (0.02-1.5 Mq at la limits). 


Calchi Novati et al.l ( 20051) present an analysis of the POINT-AGAPE data. They 
indicate that their microlensing event rate is much larger than self-lensing expectation 
alone. This leads to a lower limit of 20% of the halo mass in the form of MACHOs with 
mass between 0.5-1 Mq at the 95% conhdence level. The lower limit drops to 8% for 
MACHOs ~ 0.01 Mq. 


On the contrary, the MEGA team (jde Jong et al 


20061) conclude that their 14 events 


are consistent with self-lensing prediction. This rules out a MACH O halo fraction large r 


than 30% at the 95% conhdence level. However, recent studies by 


Ingrosso et al. 


(12006 . 


20071) show that, when compared with the timescales, maximum huxes, and spatial 


distributions from Monte-Carlo simulations, the MEGA events cannot be fully explained by 
self-lensing alone. 

The evidence for MACHOs can be inferred from individual events provided good 
sampling of the light curves. For example, the WeCA PP-GLl/POINT-AGAPE-S3 e vent 


was independe ntly discoverec 


the WeCAPP ( Riffeser et al 


by t he POINT-AGAPE ( Panlin-Henriksson et al. 


2003) and 


20031 ) collaboration. Joint analysis of the light curves from 


these two collaborat ions has led to the conclusion that this event is very unlikely to be 


a self-lensing event (IRiffeser et al. 


20081) . This is because with a realistic model of the 


three-dimensional light distribution, stellar population and extinction of M31, an self-lensing 
event with the parameters of WeCAPP-GLl/POINT-AGAPE-S3 is only expected to occur 
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every 49 years. In contrast, a halo lensing event (with 20% of the M31 halo consists of 1 
Mq MACHOs) would occur every 1 0 yea rs. In addition, combining the data from PLAN 


and WeCAPP. 


OAB-N2. 


Calchi Novati et al 


Calchi Novati et ah 


2010l) have revealed the hnite-source effect in the event 


I 2 OIO ) then used the hnite-source effect to determine the size 


of the Einstein ring radius, and combined the Einstein ring radius with the time required 
for the lens to travel the Einstein ring radius {t^) to derive the lens proper motion. Hence 
they were able to put a strong lower limit on the lens proper motion. Since the lens proper 


motion in the halo lensing scenario is different (anc 


scenario (see e.g. Fig. 4 of 


Calchi Novati et al 


much larger) than the self-lensing 


2010h . their result thus favors the MACHO 


lensing scenario over self-lensing for this event. 


To conclude, some statistical studies and individual microlensing events point to a 
non-negligible MACHO population, though the fraction in the halo mass remains uncertain. 
To pinpoint the MACHO mass fraction, a better understanding of the luminous population, 
i.e. self-lensing rate, is needed. 


Determining the MACHO mass fraction requires precise detection efficiency studies 
which is beyond the scope of this publication. We have started such efficiency studies 
(Koppenhoefer et ah, in preparation). By simulating artihcial events into the WeCAPP 
images we will be able to account for true noise and systematics introduced by the processing 
and the detection procedure. The results will enable us to derive a robust estimate of the 
MACHO mass fraction using an improved M31 model for disk, bulge and halo. 


7. Summary and Outlook 

WeCAPP has monitored the bulge of M31 for 11 years, which is the longest time-span 
among M31 microlensing surveys. We have established an automated selection pipeline and 
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present the final 12 microlensing events from WeCAPP. The brightest event of it (see e.g. 


Riffeser et al. 


2008f) is hard to reconcile with self-lensing alone, hence hints at the existence 


of MACHOs in the halo of M31. A similar case has been found by Calchi Novati et al. 

( 2010 ). 

To gain insights to the MACHO fraction, in-depth understanding of self-lensing is 
required, which is only possible with microlensing experiments that monitor both the bulge 
and the disk of M31 simultaneously. With the advent of Pan-STARRS 1, we have conducted 
a dedicated project monitoring M31 with a £eld-of-view of ~ 7 degree^. Prelini inary results 


of six microlensing events in the bulge have been published in 
analysis will follow soon (Seitz et al. in prep.). 


Lee et al 


( 20121) . and a full 
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9. Appendix 

We show the degenercy between fpwHM and Ap for our 12 WeCAPP microlensing 
events in Fig. [13 These two parameters are highly degenerate for very short time-scale 
events (tpwHM below 2 day). For events with long time-scale, i.e. W3, W4, W5, W9, WIO 
and Wll, their tpwHM and Ap are well constrained. 

We also show the degeneracy between tp and uq for our 12 WeCAPP microlensing 
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events in Fig. [161 In the pixel-lensing regime, and uq are highly degenerate, hence the 
tE-Uo maps are much more degenerate than the tpwHM-^F maps. 
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Season 

Observatory 


R-band 



I-band 


R 

or I-band 

R or I-band 



FI 

F2 

F3 

F4 

FI 

F2 

F3 

F4 

FI 

F2 

F3 

F4 

FI or F2 or F3 or F4 

1997 1998 

WS 

36 

7 

1 

4 

33 

7 

0 

3 

37 

7 

1 

4 

38 

1998 1999 

WS 

33 

1 

1 

1 

28 

1 

1 

1 

33 

1 

1 

1 

33 


WS 

64 

0 

0 

0 

60 

0 

0 

0 

65 

0 

0 

0 

65 


CA 

89 

0 

0 

0 

84 

0 

0 

0 

95 

0 

0 

0 

95 

1999 2000 

WS or CA 

128 

0 

0 

0 

124 

0 

0 

0 

134 

0 

0 

0 

134 


WS 

75 

0 

16 

0 

68 

0 

15 

0 

75 

0 

16 

0 

75 


CA 

106 107 107107 

89 

89 

89 

89 

106 107107107 

107 

2000 2001 

WS or CA 

153 107119 107 

137 

89 101 

89 

153 107119 107 

154 


WS 

106 

0 

23 

0 

93 

0 

21 

0 

106 

0 

23 

0 

106 


CA 

134 136 136 136 

119 119 119 119 

137137137137 

137 

2001 2002 

WS or CA 

200 136 148 136 

176 119 129 119 

201 137 148 137 

201 


WS 

27 

11 

16 

11 

24 

10 

17 

12 

27 

11 

18 

12 

45 


CA 

7 

7 

7 

7 

6 

6 

6 

6 

7 

7 

7 

7 

7 

2002 2003 

WS or CA 

34 

18 

23 

18 

30 

16 

23 

18 

34 

18 

25 

19 

52 

2003 2004 

WS 

35 

24 

29 

31 

33 

21 

26 

29 

35 

24 

29 

32 

69 

2004 2005 

WS 

25 

23 

26 

25 

19 

16 

19 

19 

26 

23 

26 

25 

47 

2005 2006 

WS 

30 

25 

28 

28 

26 

20 

22 

23 

32 

26 

28 

28 

71 

2006 2007 

WS 

107106 103 103 

48 

45 

46 

47 

107 108 104 103 

124 

2007 2008 

WS 

62 

56 

52 

58 

36 

35 

35 

38 

63 

58 

55 

61 

92 

total 

WS 

600 253 295 261 

468 155 202 172 

606 258 301 266 

765 

total 

CA 

336 250 250 250 

298 214 214 214 

345 251 251 251 

346 

total 

WS or CA 

843 503 530 511 

690 369 402 386 

855 509 536 517 

1015 


Table 1: The number of analyzed nights per year during the 11 WeCAPP seasons. Note that 
from 1999 until 2002 we used both telescopes at Wendelstein (WS) and Calar Alto (CA). A 
season is dehned to last from May 1*^ until 30*^ April of the next year. The total amount 
of observed nights are 1015 nights out of 11 years. A total of 4,432 stacked frames were 
analyzed in both hlters and 4 helds. 
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criterion number left LCs from I colors in Fig. [H] 


I 

Analyzed LCs with > 50 data points 

3,872,240 

100.0% 


II 

Three consecutive 3a in R 

719,628 

18.6% 


III 

Xr < 1.5 and xi < 2-1 

152,753 

3.9% 

yellow 

IV 

{S/N) > 8.9 in R and maximum with good PSF 

2,247 

0.5% 

black 

V 

-1.5 > < 1.5 

1,379 

10.4% 

green 

VI 

^FWHM < lOOOd 

63 

16.8% 

blue 

VII 

samp^^^° >0.18 and samp^^^° >0.08 and 
samp^^^° >0.18 and sampY^° >0.08 

12 

9.3% 

red 


Table 2: Selection criteria for the 11 year data taken between MJD = 685.5 and 4535.3. In 
the 3rd column we show how applying each criterion (II, III, IV, V, VI, VII) to the 3,872,240 
light curves (criterion I) reduces the hltered light curves. 


Criteria 

Number 

Percentage 

(1) I II 

719628 

100.00 

(2) I II III 


152753 

21.23=(2)/(1) 

(3) I II 

IV 

18803 

2.61=(3)/(1) 

(4) I II 

V 

402746 

55.97=(4)/(l) 

(5) I II III 

V 

106200 

14.76=(5)/(1) 

(6) I II 

IV V VI VII 

1338 

0.90=(10)/(6) 

(7) I II III 

V VI VII 

41214 

0.03=(10)/(7) 

(8) I II III IV VI VII 

15 

80.00=(10)/(8) 

(9) I II 

IV VIVII 

5382 

0.22=(10)/(9) 

(10) I II III IV V VI VII 

12 


Table 3: Number of events that pass a subset of the criteria. 













RA(J2000) 

Dec(J2000) 

Am31 

to 

Rwhm 

rriR 

^fr 

mj 

Api 

(R-I) 

log(i£;) 

log(wo) 

Xlof 




[arcmin] 

MJD [day] 

[day] 

[mag] 

[10-5 Jy] 

[mag] 

[10-5 Jy] 

[mag] 

[day] 

[^e] 


WeCAPP-Olf 

00:42:30.03 

41:13:01.5 

4.08 

1850.84^°;“ 


18.68 

10.29t°:54 

17.81 

1 S OO+0-99 
lO.zz_o.98 

0.88l[;;°5 

7 oc+0.01 

' -o^-O.Ol 

Q 1 O+0.02 

— O.lO_ 0,02 

1.21 

WeCAPP-02 

00:42:33.01 

41:19:58.5 

4.36 

1895.4lt°:7i 

5.85t^:°^ 

20.81 

i.45t°:57 

19.64 

Q 00+0.88 

O.OO_ 0.66 

1 1 T'+O.OO 
-‘-•'*-'-0.09 

1 61+°-29 

i.OO- 0.49 

1.20 

WeCAPP-03 

00:42:57.03 

41:12:27.9 

4.40 

1585.35;[][8 

^•^^-0.26 

20.47 

1 qq+O-iO 
-‘-•^^-0.04 

19.11 

rc Cl +0.26 
^•^-^-0.11 

1 4A+0-04 

4-00_o.o4 

2-641^62 

—2 7q+^-®^ 
' ^-0.62 

1.23 

WeCAPP-04 

00:42:54.07 

41:14:37.0 

2.48 

2178.86;^^ 

'^•^^-0.27 

20.18 

^•^^-0.10 

19.58 

'^•^’^-0.16 

0-60l°;95 

0.99l°;ll 


1.40 

WeCAPP-05 

00:43:02.03 

41:18:29.2 

4.09 

2177.85;°j9 

a 04+1.06 

20.94 

-‘-•^^-0.09 

20.38 

1 71 +0.20 
' -^-0.16 

0.561^11 

1 47+0.33 
-0.32 

1 1 q+0.39 
J^-^^-0.42 

1.38 

WeCAPP-06 

00:42:49.01 

41:14:55.3 

1.52 

2317.23t[];!]| 

0 44+0.16 

18.15 

16.794;^^ 

17.36 

27.50l^j^ 

n 7q+0-03 
^^-0.04 

Q 9'7+0.02 
-0.02 

-9.181^4 

1.67 

WeCAPP-07 

00:42:55.03 

41:18:50.9 

3.38 

1847.44;°;)^ 


20.47 

1 qq+0'3'^ 
-L-yy_o.30 

19.74 

3.08t°;59 

n 74+0.13 
^ ^-0.13 

8 . 681 °;°? 

_q 14+0.05 
^•'‘^^-0.06 

1.02 

WeCAPP-08 

00:42:50.03 

41:18:40.6 

2.75 

2111.55t°:J° 

o.584]:[6 

19.11 

6-93t[S:«I 

18.36 

10.9711:59 

o.75l[]:!]5 

i-oil?;^! 

1 ^^+0.39 
'-0.38 

1.17 

WeCAPP-09 

00:42:44.07 

41:12:54.9 

3.22 

2231.67t[!j? 

14.04t)J| 

21.17 


20.48 

1 cc+0.11 
-‘-•^^-0.10 

0.6819;!]! 

i ro+0.14 
-^■^'^-0.15 

-0.86l°;9° 

1.42 

WeCAPP-10 

00:42:12.01 

41:09:21.5 

9.07 

4026.88;[][i^ 

9 7c:+0-42 
' ^-0.48 

21.01 

1 20“^*^'^^ 
09 

20.10 

0 91 +0.26 
^•^'‘^-0.22 

0.92l°;ll 

0.811°;°° 

-0-83l°:°° 

1.38 

WeCAPP-11 

00:42:46.02 

41:15:11.9 

1.01 

1198.8ltS];[^ 

6-4i4S;^°9 

21.09 

1 1 0 +O.O 8 

i.4O_0.07 

19.89 

2.691^16 


1 91 +0.47 
-^■^-'--0.36 

- 0 . 88 l°;« 

1.62 

WeCAPP-12 

00:43:07.08 

41:17:40.7 

4.64 

4018.87;o;2i 

4.92t[;°8 

20.64 

1.69t°:50 

19.80 

2.901!;;55 

o-84l9:12 

i.UZ_o.23 

-0.79l!;;ll 

1.17 


Table 4: Paczynski parameters of the WeCAPP microlensing events. In the fourth column we also show the distance 
from the M31 center. ^To determine the parameters of the light curve, we include data from POINT-AGAPE, as in 
Riffeser et ah (2003). 
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Fig. 10.— WeCAPP microlensing event light curves: WeCAPP 1-4 with corresponding cut¬ 
outs of the difference frames in R and I band. The data points are color-coded in grey-scale 
according to their errors; measurements with larger errors in R (I) are shown in light-red 
(light-blue), while measurements with smaller errors are shown in dark-red (dark-blue). 
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Fig. 11.— WeCAPP microlensing event light curves: WeCAPP 5-8 with corresponding 
cut-outs of the difference frames in R and I band. 
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Fig. 12.— WeCAPP microlensing event light curves: WeCAPP 9-12 with corresponding 
cut-outs of the difference frames in R and I band. 
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Table 5: Overview of the 56 microlensing events in M31 presented in different papers. In 
the 3rd column for the total number we counted re-published detection only once. For 
comparison see Figure [HI 
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Fig. 13.— to, tpwHM and distribution of all microlensing events in M31 reported up to 
now. Short events seem to be frequent. Around 20 days there is a gap which separates short 
from long lensing events. The WeCAPP events are plotted in red, VATT/Columbia from 
Crotts & Tomaney (1996) in violet, POINT-AGAPE in blue, Belokurov-AGAPE in cyan, 
MEGA in green, AGAPE Z1 in gray, NMS in magenta, PLAN in yellow, SLOTT-AGAPE in 
brown, VATT/Columbia from Uglesich et al. (2004) in black, PAndromeda in orange. Note 
that the above distributions can not be directly taken to study the nature of events (halo vs. 
self-lensing) and to constrain the halo MACHO fraction, but instead have to be corrected 
for the surveys’ detection efficiencies. 
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Fig. 14.— 56 microlensing events in M31. The image shows a wide held (50 x 70 arcmin) 
from the 2”'^ Palomar Sky Survey (ESO Online Digitized Sky Survey - DSS-2-blue). The 
WeCAPP events are plotted in red, POINT-AGAPE events in blue, Belokurov events in 
cyan, MEGA events in green, AGAPE Z1 event in black, NMS event in magenta, PLAN 
events in yellow, PAndromeda events in orange, Slott-AGAPE in brown, VATT from Grotts 
& Tomaney (1996) in violet, and VATT from Uglesich et al. (2004) in black. We note that 
although our event W8 appears to be on top of the OAB-N2 event in this hgure, their tg differ 
by more than 6.5 years (W8 has a to of JD=2452112 and OAB-N2 has a to of JD=2454467). 
In addition, from our difference images their positions differ by 0.5 arcsec. 
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Fig. 15.— tpwHM vs. Apji distribution of the 12 WeCAPP microlensing events. We note 
that the tpwHM and Air^-scale in the hgures are different for each event. The 1, 2, and 3-cr 
contours are marked in green, cyan and blue, respectively. The best-£t parameter is marked 


in magenta. 
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Fig. 16.— VS. Mo distribution in log scale of the 12 WeCAPP microlensing events. We 
note that the t-^ and Mo"Scale in the hgures are different for each event. The 1, 2, and 3-cr 
contours are marked in green, cyan and blue, respectively. The best-fit parameter is marked 
in magenta. 
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